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Abstract 

In  a  series  of  earlier  studies,  we  have  shown  that  a  variety  of  solution,  liquid  and 
solid  phase  properties  can  be  represented  analytically  in  terms  of  quantities  related  to  the 
electrostatic  potentials  on  molecular  surfaces.  These  quantities  include  the  positive  and 
negative  extrema,  the  positive  and  negative  average  values  and  variances,  and  the  average 
deviation.  We  have  now  investigated  how  sensitive  these  quantities  are  to  the  molecular 
conformation.  Surface  potentials  were  computed  at  the  HF/6-31G*  level  for  a  total  of  35 
conformations  of  ten  different  molecules.  Our  overall  conclusion  is  that  conformational 
effects  upon  applications  of  molecular  surface  electrostatic  potentials  are  most  likely  to  be 
of  concern  if  (a)  formation  of  the  conformer  considerably  diminishes  internal  polarity , 
and/or  (b)  the  application  in  question  is  strongly  dependent  upon  the  positive  variance. 
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statistical  quantities 


Introduction 


The  electrostatic  potential  V(r)  created  in  the  space  around  a  molecule  by  its  nuclei 
and  electrons  is  well  established  as  a  guide  to  molecular  reactive  behaviour.  (For  reviews, 
see  refs.  [1-6].)  It  is  defined  by  eq.  (1),  in  which  the  molecule  is  treated  as  a  collection  of 
stationary  point  charges,  the  nuclei,  surrounded  by  a  continuous  but  static  distribution  of 
electrons: 


V(r)  =  X 

A 


Za 

RA-r| 


Cp(r')dr' 

J  |r'  -r| 


is  the  chtirge  on  nucleus  A,  located  at  R^,  and  p(r)  is  the  electronic  density  function  of 

the  molecule.  The  first  term  on  the  right  side  of  eq.  (1)  is  the  nuclear  contribution  to  V(r), 
and  is  positive;  the  second  term  is  due  to  the  electrons  and  is  accordingly  negative. 

Over  a  period  of  many  years,  the  electrostatic  potential  was  used  extensively  as  a 
tool  for  identifying  and  ranking  the  molecular  regions  most  susceptible  to  electrophilic 
and/or  nucleophilic  attack  and  for  determining  general  patterns  of  positive  and  negative 
potential  that  promote  or  inhibit  molecular  interactions,  such  as  those  between  drugs  and 
receptors.  Bernard  and  Alberte  Pullman  and  their  collaborators  were  pioneers  in  applying 
the  electrostatic  potential  to  the  analysis  of  particularly  biochemical  systems;  some  of  their 
numerous  contributions  are  discussed  in  refs.  [2-4].  The  quantitative  analysis  of  V(r) 

initially  emphasized  locating  and  evaluating  the  most  negative  potentials,  These  are 

usually  associated  with  (a)  the  more  electronegative  atoms,  such  as  N,  O,  F,  Cl,  S  and  Br, 
and  (b)  unsaturated,  aromatic  and  strained  carbon-carbon  bonds.  The  magnitudes  of  the 
Vjnin  can  often  be  related  to  reactive  properties,  for  instance  the  pK,  values  of  azine 
nitrogens  [5],  and  epoxide  carcinogenicity  [4].  More  recently,  attention  has  focused  upon 
the  electrostatic  potential  Vs(r)  on  the  molecular  surface.  It  was  shown  that  both  the 

positive  and  negative  surface  extrema.  Vs  max  ^s,min’  related  to  reactive 
behavior.  For  example,  Vg^max  and  Vg  min  for  a  large  variety  of  molecules  correlate  with 
hydrogen  bond  acidity  and  basicity,  respectively  [7].  (We  take  the  molecular  surface  to  be 
the  0.001  au  contour  of  the  electronic  density,  as  suggested  by  Bader  et  al  [8].  We  have 

earlier  discussed  the  use  of  p(r)  contours  to  define  molecular  surfaces  [9,10].)  This 
approach  was  still  limited  in  scope  however;  Vmin.  Vg  min  and  Vg  max  are  certainly  key 
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features  of  the  molecular  electrostatic  potential,  but  these  site-specific  quantities  do  not 
convey  all  the  information  that  is  contained  in  V(r). 

Accordingly,  in  recent  years  we  have  sought  to  develop  mechanisms  for  more 
adequately  describing  and  quantitatively  characterizing  the  electrostatic  potential  over  an 
entire  molecular  surface.  We  have  found  that  this  can  be  achieved  through  the  introduction 
of  several  statistically-defined  global  quantities  that  explicitly  reflect  the  magnitude  of  VgCr) 
at  each  point  on  the  molecular  surface  [11-14]: 

(a)  Vg  ,  Vg  and  Vg  are  the  positive,  negative  and  overall  average  potentials  on  the 


surface. 


v^=-ivs>,) 

—  1  ^ 

Vs  - - SVf(r,) 

“1=1 


1  r  ^  m 

Vs=—  IVsV.)  +  XVs'Cr.) 


m-l-n 


(b)  n  is  the  average  deviation  of  Vs(r), 


,  m+n  _ 

n - £|Vs(ri)-Vs| 

m  -L  n  '  ' 


m+n  ' 

which  we  interpret  as  a  measure  of  the  local  polarity,  or  internal  charge  separation,  that 
is  present  even  in  molecules  with  zero  dipole  moment. 

(c)  a+,  ai  and  are  the  positive,  negative  and  total  variances  of  Vgfr),  which  reflect 


the  range  or  variability  of  VgCr),  emphasizing  its  extrema. 


<^i  =  -i[vsVi)-Vs^] 

(6) 

^  2 

1  m  ^ 

=  Vs(r.)-Vs- 

(7) 

2  2  2 

(8) 

(d)  V  indicates  the  degree  of  balance  between  the  positive  and  negative  surface  potentials. 
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V  = 


g+g- 


(9) 


When  g+  =  gi,  then  v  achieves  its  maximum  possible  value  of  0.25. 

In  a  series  of  studies,  we  have  demonstrated  that  it  is  possible  to  develop 
quantitative  analytical  relationships  of  good  accuracy  for  a  variety  of  solution,  liquid  and 

solid  phase  properties  in  terms  of  some  subset  of  these  quantities  (plus  ^inin>  ^S.max  ^tld 

Vs  niin),  which  are  computed  for  the  individual,  isolated  molecules.  These  properties 
include  heats  of  fusion  [15],  vaporization  [11]  and  sublimation  [16],  solubilities  [17-19], 
solvation  free  energies  [20,21],  boiling  points  and  critical  constants  [22],  partition 
coefficients  [23,24],  diffusion  coefficients  [25],  surface  tensions  [15],  lattice  energies  [26], 
liquid  and  solid  densities  [  1 5] ,  and  impact  sensitivities  [27] . 

In  calculating  the  quantities  defined  by  eqs.  (l)-(9),  our  first  step  is  to  optimize  the 
molecular  geometry,  so  that  we  are  presumably  dealing  with  the  most  stable  conformation. 
An  important  and  frequently-posed  question  concerns  the  sensitivity  of  the  computed 
quantities  (and  subsequent  relationships)  to  the  conformation  of  the  molecule.  How  much 
is  the  surface  electrostatic  potential  affected  if  the  molecule  is  induced  to  adopt  another 
conformation,  perhaps  due  to  thermal  or  environmental  factors?  Our  objective  in  the 
present  work  has  been  to  address  this  question. 


Procedure  and  Results 

For  the  molecules  1-10,  we  have  investigated  the  conformations  shown  in  Figures 
1  and  2.  Each  molecular  geometry  was  first  fully  optimized  in  the  ground  state,  and  then 
re-optimized  in  one  or  more  different  conformations  (local  minima).  All  calculations  were 
carried  out  with  Gaussian  92  [28]  at  the  HF/6-31G*  level.  For  each  conformation,  we 

computed  the  surface  area,  Vs,max  and  Vs^^in.  and  the  properties  defined  by  eqs.  (2)-(9). 

The  results  are  listed  in  Table  1 ,  along  with  the  HF/6-3 1 G*  relative  energies  within  each 
conformer  group. 

CH3-CH3  CH3-CH2F  CH3-CH2NO2  CH3-CH3OH  CH3CI-CH3F 

1  2  3  4  5 

CH2NH2-CH2NO2  CH2OH-CH2OH  CH3-CH2-CH2OH  CH2=CH-C00H  NH=CH-CH2N0 

6  7  8  9  10 
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Discussion 

Ethane  and  Its  Monosubstituted  Dervatives 

Both  the  staggered  and  the  eclipsed  conformer  have  been  examined  for  molecules 
1-4.  The  latter  is  invariably  the  less  stable,  by  approximately  3  -  4  kcal/mole.  (The 
experimental  value  for  ethane  is  3.0  kcal/mole  [29].)  The  computed  properties  tend  to  be 
quite  similar  for  each  pair  of  conformers.  In  the  case  of  4,  we  also  considered  the 
structure  (4c)  resulting  from  a  64°  rotation  around  the  C-0  bond  of  the  staggered  form. 
This  requires  an  energy  of  only  0.1  kcal/mole,  but  produces  a  rather  large  change  in  the 

global  property  V^,  indicating  its  sensitivity  to  the  precise  juxtaposition  of  the  hydroxyl 
and  neighboring  hydrogens. 

1, 2-Disubstituted  Ethanes 

Four  conformers  were  investigated  for  each  disubstituted  ethane.  The  least  stable 
structure  is  usually  that  in  which  the  two  substituents  are  eclipsed;  8c  is  an  exception, 
presumably  due  to  the  short  distance  (1.95  A)  between  one  of  the  methyl  hydrogens  and 
one  on  the  substituted  carbon.  In  terms  of  our  present  objective,  it  is  5d  that  is  of 
particular  interest,  in  that  it  differs  significantly  from  its  conformers  in  terms  of  nearly  all 
of  its  computed  surface  properties.  This  is  not  observed  for  any  other  of  the  disubstituted 
ethanes.  In  5a,  5b  and  5c,  the  negative  halogen  potentials  evidently  overlap  and  reinforce 
each  other,  creating  a  relatively  strong  negative  region  and  a  corresponding  positive  one. 

In  5d,  on  the  other  hand,  the  negative  region  are  separate  and  therefore  weaker.  The 
magnitudes  of  n  confirm  that  5d  has  much  less  internal  polarity  than  do  its  conformers. 
For  6-8,  the  surface  properties  generally  tend  to  be  fairly  similar  among  the  conformers. 

The  most  striking  exceptions  to  this  are  a+,  a_  for  8a  and  Vs  for  8b. 

Molecules  9  and  10 

The  four  conformers  of  9  fall  neatly  into  two  groups;  9a  and  9c  are  significantly 
less  stable  than  9b  and  9d,  presumably  because  of  the  proximity  of  the  hydrogens  in  the 
former,  which  creates  considerable  internal  polarity.  The  surface  properities  are  fairly 
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uniform  within  the  two  groups.  Six  different  conformers  of  10  have  been  investigated. 
Overall,  they  do  not  vary  markedly  in  their  surface  properties. 

Summary 


From  the  data  in  Table  1,  certain  generalizations  can  be  made  concerning  the  effects 
of  conformational  changes  upon  the  molecular  surface  properties  of  present  interest: 

(a)  The  surface  areas  are  only  slightly  affected. 

(b)  The  two  site-specific  properties,  Vg  and  Vg  and  the  global  properties  Vg 

and  n  usually  change  relatively  little.  The  only  exceptions  to  this  are  5d  and  the 

pair  9a/9c,  which  were  formed  by  rotations  that  either  eliminated  or  produced  regions 
of  markedly  reinforced  negative  or  positive  potential.  This  is  then  reflected  in  all  of  the 

properties  of  Y^{r).  However  increasing  the  internal  polarity  can  normally  be  expected 

to  significantly  decrease  molecular  stability,  so  that  conformers  such  as  9a  and  9c 
are  less  apt  to  play  important  roles. 

2 

(c)  G+  is  the  most  sensitive  to  conformational  variations.  It  is  probable  that  V  andv 
will  consequently  be  affected,  although  perhaps  not  to  the  same  extent  because 
they  also  include  a_. 

On  the  basis  of  the  results  in  Table  1,‘ therefore,  it  appears  that  conformational  effects  upon 
applications  of  molecular  surface  electrostatic  potentials  are  most  likely  to  be  of  concern  if 
(a)  formation  of  the  conformer  considerably  diminishes  internal  polarity,  and/or  (b)  the 

application  in  question  is  strongly  dependent  upon  c + . 
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Figure  Captions 

Figure  1.  Conformers  investigated  for  molecules  1  -  6. 
Figure  2.  Conformers  investigated  for  molecules  7  - 10. 
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